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We report the measurement of an inclusive partial branching fraction for charmless semileptonic 
B decay and the extraction of \V u b\- Candidates for B — > X u iv are identified with a novel X u - 
reconstruction method based on neutrino reconstruction via missing 4-momentum and a technique 
called "simulated annealing." Based on 86.9 fb -1 of data taken with the Belle detector, we obtain 
AB(B -» X u £p;Mx < 1.7 GeV/c 2 ,g 2 > 8.0 GeV 2 /c 2 ) = (7.37 ± 0.89(stat.) ±1.12(syst.) ±0.55(6 -> 
c)±0.24(6 -» u))xl0 -4 and determine \V ub \ = (4.66±0.28(stat.)±0.35(syst.)±0.17(6 -> c)±0.08(6 -> 
u) ± 0.58(theory)) x 10" 3 . 

PACS numbers: 12.15.Hh, 13.25.Hw, 29.85.+C 



The off-diagonal element V u b in the CKM matrix plays 
an important role in CP-violation and rare decays of the 
B meson. It is an important ingredient in overconstrain- 
ing the unitarity triangle by measuring its sides and an- 
gles. In the experiments on the T(4S f ) resonance, its mag- 
nitude is extracted from measurements of the B — > X u lv 
process in the limited region of lepton momentum [l| 
or the hadronic recoil mass Mx [2| where the contribu- 
tion of background from the B — > X c tv process is sup- 
pressed. These experiments achieve more precise mea- 
surements than LEP experiments [3( due to higher signal 
purity; however, the need to extrapolate measured rates 
from such limited regions results in large theoretical un- 
certainties on \V u b\- A recent theoretical development 
suggests that one can significantly reduce the theoretical 
uncertainty on the extrapolation by applying simultane- 
ous cuts on Mx and the invariant mass squared of the 
lepton-neutrino system (q 2 ) in inclusive B — > X u lv Q. 
We report here the first result with simultaneous require- 
ments on Mx and q 2 . The result is obtained with a 
novel X„-reconstruction method based on a combination 
of neutrino reconstruction and a technique called simu- 
lated annealing |5j to separate the two B meson decays. 
This method allows us to measure Mx and q 2 with good 
efficiency so that it achieves good statistical precision and 
small theoretical uncertainty with a modest integrated 
luminosity. This analysis is based on 78.1 fb _1 data, cor- 
responding to 85 million BB pairs, taken at the T(4S f ) 
resonance, and 8.8 fb" 1 taken at an energy 60 MeV be- 
low the resonance, by the Belle detector Iff at the energy- 
asymmetric e + e~ collider KEKB Q. 

We select hadronic events containing one lepton can- 
didate (electron or muon) having momentum above 



1.2 GeV/c in the center-of-mass-system (CMS) of the 
T(4S f ). To remove events with more than one neutrino, 
we exclude events containing additional lepton candi- 
dates (p* > 0.5 GeV/c and p* > 0.8 GeV/c). The neu- 
trino is reconstructed from the missing 4-momentum in 
the event (p„ = p T (4S) - ^iPi, E u = ^r(4S) ~ 
The net observed momentum (Ejpi) and energy (Xj_E!j) 
are calculated using particles surviving track quality cuts 
based on the impact parameter to the interaction point 
and a shower energy cut. Pairs of pions and electrons 
passing secondary vertex criteria are treated as Kg and 
photons, respectively. All remaining charged tracks are 
classified as kaons, pions, or protons, based on parti- 
cle identification information. K\ candidates are iden- 
tified from isolated clusters of hits in the detector for 
K^s and muons [8j. The energy of each particle can- 
didate is calculated based on its momentum and mass 
assignment. We then compute the missing mass of 
the event, defined as MM 2 = E 2 /^ — \p„\ 2 /c 2 , where 
the sign of E 2 is reversed when E v < 0. We require 
-1.5 GeV 2 /c 4 < MM 2 < 1.5 GeV 2 /c 4 to suppress events 
with missing particles and with particles removed due to 
poor reconstruction quality. For events that pass this 
requirement, we add back tracks and clusters rejected 
earlier due to reconstruction quality, selecting the com- 
bination that gives the smallest value of |MM 2 |. This 
determines the set of particles that are used in the sub- 
sequent analysis. Events are further required to have a 
net charge of or ±1, and a polar angle for the missing 
momentum within the barrel region (32° < 8 < 128°). 
To suppress beam-gas events, we demand that the net 
charge of all proton candidates be zero. Requiring that 
the cosine of the angle of E\ candidates with respect to 
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FIG. 1: Distributions for 3 discriminants and W before 
(dashed histogram) and after (solid histogram) simulated an- 
nealing, for real data. Distributions for 3 discriminants for 
correct combination of particles for MC events (solid curve). 



the missing momentum be less than 0.8 rejects events 
where the neutrino candidate is actually a K L meson. 

We then seek the most likely combination of parti- 
cles belonging to Xlv, the remainder being from the 
associated B -meson (_Z? pp)- Six discriminant variables 
are used: the momentum, energy, and polar angle of 
B opP (p*b, E* Bl cos9* B ) in the CMS, its charge multi- 
plicity (-/V c h) and net charge times the lepton charge 
{Qb x Qe), and the missing- mass squared recalculated 
with the energy and mass of -B pp constrained to the 
known values (MM\ lv ). Using Monte Carlo (MC) sim- 
ulation events for T(45) — > BB where at least one B 
decays into Xlv, we determine probability density func- 
tions (PDFs) for correct Xiv combinations and for ran- 
dom Xlv combinations. Random candidates for Xlv 
consist of the lepton and neutrino candidates plus parti- 
cles from the remainder of the event, selected randomly 
so that the relative multiplicities between X and the re- 
mainder of the event matches that at the generator level. 
From the PDFs we calculate two likelihoods, £ (correct) 
and £ (random). The most likely candidate combina- 
tion in each event is found by minimizing the parameter 
W = £(random)/(£(random) + ^(correct)). 

To minimize W, we have developed an approximate 
iterative algorithm based on simulated annealing. We 
start from the initial candidate for Xlv that consists of 
the lepton and neutrino plus approximately one third 
of the remaining particles, selected randomly. We move 
a random particle (other than the lepton or neutrino) 
between the X and i? pp sides in an iterative way, where 
in one iteration we cross all particles at least once, and 
search for the combination that gives the minimum W 
with 50 iterations. During the iteration process we take 
special care to reduce the chance of convergence to a local 
minimum of W . For instance, after every fifth iteration 
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FIG. 2: Mx (a) and q 2 (b) resolution distributions for 
B — > X u lv MC events. Histograms (curves) show the results 
with the simulated annealing method (with correct particle 
assignment to X u ). 
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FIG. 3: (a) M x distributions for B 
■ 2 distributions for B 



D*lv control samples, 
(b) q z distributions for B — » J/tpX control samples. Points 
are data and the histogram is MC. 



we compare all combinations that can be constructed by 
crossing one particle and use the combination that gives 
the largest value of W to seed a new cycle. We repeat 
this iteration process 10 times, starting each time with 
a different initial candidate, and select the case with the 
smallest W. Figure shows the distributions in three of 
the six discriminant variables and W, before and after 
simulated annealing. Also shown are the distributions 
for the correct combination in signal MC events. 

The final candidate is required to satisfy: i) W < 0.1, 
ii) 5.1 < E* B < 5.4 GeV, iii) 0.25 < p* B < 0.42 GeV/c, 



-0.2 < MM 



Xlv 



< 



iv) -2 < x Q B < +1, and 
0.4 GeV 2 /c 4 . Contamination from the continuum is re- 
duced by demanding \ cos 9 bi\ < 0.8, where 9 Be is the 
angle between the thrust axis of i? pp and the lepton 
momentum. Figure [21 shows the resolutions in Mx and 
q 2 for B — > X u lv MC events. Also shown in the figure 
are the resolutions for correct combination of particles. 

The validity of the method is checked with two data 
samples, one containing 38,600 fully-reconstructed B — » 
D*iv decays and the other containing 84,100 B — > 
J/ipX, J /if) — > l + l~ decays. These data samples are also 
used to calibrate the detection efficiency. For the J/i/jX 
sample we treat one of the two leptons from J ftp as a neu- 
trino, to emulate B — > Xiv. Corresponding MC events 
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FIG. 4: (a) M x distribution for q 2 > 8.0 GeV 2 /c 2 . (b) q 2 
distribution for Mx < 1-7 GeV/c 2 . Points are the data and 
histograms are backgrounds from D*tv (dotted), Dlv (short 
dashed), others (long dashed), and total background contri- 
bution (solid). Lower plots show the data after background 
subtraction. Solid curves show the inclusive MC predictions 
for B -> X u lu. 



are generated with the QQ event generator [9| a nd the 
detector response is simulated using Geant 3 10] . Fig- 
ure[3]shows the Mx distribution for the D*£ v sample and 
the q 2 distribution for the J/ i/jX sample, with MC distri- 
butions scaled to the number of background-subtracted 
events in the respective data samples. The peaks are at 
the D* mass and J/ ip mass squared, as expected, and the 
shapes are in good agreement between data and MC. Al- 
though these results verify that the simulated annealing 
method works as expected, we observe a small difference 
of efficiency between the data and MC. By averaging the 
results of the two data samples, we obtain the efficiency 
ratio r cff = 0.891 ± 0.043 between the data and MC. 

We observe 8910 events in the X u lv "signal" region, 
defined as M x < 1.7GeV/c 2 and q 2 > 8.0GeV 2 /c 2 . 
These consist of semileptonic decays, B — » X CtU £v, other 
BB background, and residual continuum events. The 
continuum contribution is estimated from off-resonance 
data to be 251 ± 48 events and is subtracted directly 
from the analyzed distributions. The contributions from 
B — > X c £v and other BB backgrounds are estimated 
via MC in the "background" region M x > 1.8 GeV/c 2 , 
where X c lv dominates, and extrapolated to the signal 
region. We estimate them by fitting the Mx and q 2 
distributions from MC events to those from the data 
using a two-dimensional \ 2 fit method. Contributions 
to X c lv come from D^Hv, D**lv, and D^iriv. Their 
branching fractions are floated in the fit. The total rate 
for other BB backgrounds, arising from sources such as 
b —> c — > s£v and fake leptons, which amounts to less 
than 1% of events in the signal region, is floated. The 
small X u lv contribution is estimated iteratively and is 
found to be (0.94 ± 0.04)% of events in the background 
region. The obtained branching fractions for the exclu- 



sive B —* X c lv modes are consistent with the PDG val- 
ues 01 • The BB background in the signal region is 
estimated to be 7283 ± 130 ± 63 events, where the first 
and second errors come from fit and MC statistics, re- 
spectively. The upper plots in Figure 0] show the Mx 
distribution for q 2 > 8.0 GeV 2 /c 2 and the q 2 distribu- 
tion for Mx < 1.7 GeV/c 2 after continuum subtraction. 
After subtracting the BB backgrounds, we obtain dis- 
tributions for the X u lv signal, shown in the lower plots. 
The net signal is N ^ s = 1376±167 events where the error 
is statistical only. 

In order to extract the partial branching fraction AB 
for B — > X u lv in the signal region, a Monte Carlo simula- 
tion is used to convert iV bs to the true number of signal 
events produced in this region, iVt ruo , and to estimate 
the efficiency for these events to be observed anywhere, 
Esignai- In the MC simulation, B — > X u tv decays are 
simulated based on the prescription of ref. That 
analytic result gives 0(a s ) corrections to leading order 
in the heavy-quark expansion for the triple differential 
B — > X u lv rate and includes the effect of the 6-quark's 
Fermi motion. Two parameters therein, the 6-quark pole 
mass, mb, and the average momentum squared of the 
6-quark inside the B meson, /i 2 ,, are derived from the 
CLEO measurements of the hadronic mass moments in 
inclusive D — > X c lv and photon energy spectrum in B — 
X sl 0. We use m b = 4.80 ± 0.12GeV/c 2 and \i\ = 
0.30 ± 0.11GcV 2 /c 2 , which differs from CLEO's evalu- 
ation in that terms proportional to 1/mS and a 2 have 
been removed from the relation between the measured 
observables and m& and /i 2 . The MC events are gener- 
ated with the EvtGen generator ^J. X trU c is estimated 
by N tIue = N ohs xF(F = 1+N 2 /N 1 -N 3 /N 1 ). Here JVi is 
the number of events observed in the signal region and 
(N 3 ) is the number of events generated inside (outside) 
the signal region and observed outside (inside) the signal 
region. We find F = 0.938, and thus iVt ruo = 1291 ± 157. 
The efficiency e s i gna i is predicted to be 0.578%. We deter- 
mine AB by 0.5 x iV t rue/(e s ignal x r cff )/(2N B ), where r off 
is the efficiency correction factor described earlier, Nb is 
the number of BB events and the factor 0.5 is needed to 
take into account the electron and muon data: 

AB = (7.37 ± 0.89 ± 1.12 ± 0.55 ± 0.24) x 10~ 4 . 

The errors are statistical, systematic, from B — > X c iv 
model dependence, and B — » X u iv model dependence, 
respectively. Sources of systematic uncertainty include 
signal MC statistics (1.8%), lepton identification (2.6%), 
uncertainty of F due to imperfect detector simulation 
(1.2%), selection and reconstruction efficiency (4.9%), 
BB background estimation (14.0%). The uncertainty 
in the BB background estimation is from MC statistics 
(4.6%) and distortion of the Mx and q 2 distributions due 
to imperfect detector simulation (13.2%). Major con- 
tributions to the last error are from contamination 
(8.6%), electromagnetic cluster finding efficiency (8.2%), 
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lepton efficiency (3.2%), clusters produced by charged 
tracks (2.9%), lepton fake rate (2.7%), and K/ir sepa- 
ration (2.7%). The error from contamination is esti- 
mated using inclusive Kg events where we discard KgS to 
emulate inclusive if£ events. The error from the cluster 
finding efficiency is estimated by reducing the photon- 
finding efficiency within its uncertainty. The model de- 
pendence of X c iv is estimated to be 7.4% by varying the 
D\lv plus D\lv fraction in the D**lv by 25% and by 
varying the slope parameters of the form factors for Div 
and B*lv, p 2 D = 1.19 ± 0.19 and p 2 = 1.51 ± 0.13 [ll|, 
within their errors. B — > X u lv model dependence (3.4%) 
is estimated by varying the parameters of the inclusive 
model within their errors and by comparing to a simula- 
tion with a full exclusive implementation of the ISGW2 
model [3. 

In the context of HQET and OPE the partial branch- 
ing fraction AB(B -> X u lv) is related to \ V uh \ HE 03, 



We wish to thank the KEKB accelerator group for the 
excellent operation of the KEKB accelerator. We ac- 
knowledge support from the Ministry of Education, Cul- 
ture, Sports, Science, and Technology of Japan and the 
Japan Society for the Promotion of Science; the Aus- 
tralian Research Council and the Australian Department 
of Education, Science and Training; the National Sci- 
ence Foundation of China under contract No. 10175071; 
the Department of Science and Technology of India; the 
BK21 program of the Ministry of Education of Korea 
and the CHEP SRC program of the Korea Science and 
Engineering Foundation; the Polish State Committee for 
Scientific Research under contract No. 2P03B 01324; the 
Ministry of Science and Technology of the Russian Fed- 
eration; the Ministry of Education, Science and Sport of 
the Republic of Slovenia; the National Science Council 
and the Ministry of Education of Taiwan; and the U.S. 
Department of Energy. 



\V ub \ = 0.00444 



AB(B -> X u lv) 1.55ps 
0.002 x 1.21G(g c 2 ut ,m cut ) t b 



1/2 



cut' ' 



/hqet X ( 



4.7GoV/c s 



a 2 

ycut 



/hqet 
and 



where 1.21G(g, 

represents the fraction of events with q 2 > 
Mx < m cu t, and m^ s is one-half of the perturbative 
contribution to the mass of the T(1S). G(q 2 ut ,m cut ) 
is calculated to 0(a 2 ) and 0{\/m 2 ) in Q, including 
the effect of the Fermi motion of the b quark, which 
is expressed in terms of . We use ml = 4.70 ± 
0.12GeV/c 2 HE3, which gives G(q 2 ut ,m cut ) = 0.268 
0, 0|. The theoretical uncertainty on \V u b\ is deter- 
mined only by the uncertainty on G(q 2 ut , m cu t)- The un- 
certainty on G(q 2 ut , m cut ), in total 25%, consists of 6% for 
perturbative, 8% for nonperturbative terms (dominated 
by the weak annihilation contribution), and 23% from 
the uncertainty on m£ s 0,^1- The 23% error contains 
10% for /hqet and 13% for (mj s ) 5 . These uncertain- 
ties are positively correlated, so we add them linearly, 
whereas they have been given separately in conventional 
analyses. Using tb = 1.604 ± 0.012 ps [llj, we obtain 



\V ub \ = (4.66 ± 0.28 ± 0.35 ± 0.17 ± 0.08 ± 0.. 



x 10 



-3 



where the errors are statistical, systematic, b — > c model 
dependence, b — > u model dependence, and theoretical 
uncertainty for OPE, respectively. 

To summarize, we have performed the first measure- 
ment of I V u b I with simultaneous requirements on Mx and 
q 2 using a novel ^-reconstruction method. The result 
of I Kb I = (4.66 ±0.76) x 10~ 3 is consistent with the pre- 
vious inclusive measurements 0, 0, and the total error 
is comparable with those of the previous measurements 
on T(4S) Q, 0- Due to simultaneous requirements on 
Mx and q 2 , the /hqet error is much smaller than those 
of the previous measurements on T(45) 0,0- 
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